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In the weak disordered regime we provide analytical expressions for the electron localization 
lengths in quasi-one dimensional (Q1D) disordered quantum wire with hard wall and periodic bound- 
ary conditions. They are exact up to order W 2 (W being the disorder strength) for an arbitrary 
number of channels. Detailed numerical analysis of the Anderson localization, based on Kubo's 
formula for conductivity, show excellent agreement with analytical calculations. We establish rela- 
tionship between various lengths in Q1D systems. 
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Introduction.— A quasi-one dimensional (Q1D) geom- 
etry, as a model for a disordered wire, is of great inter- 
est in condensed matter theory. The electronic transport 
problem in weakly disordered Q1D systems can be solved 
analytically within some approximations (see, e.g., [l| for 
details). The Dorokhov-Mello-Pereyra-Kumar (DMPK) 
equation Q and random matrix theory for the transfer 
matrix Q are the two successful approaches which are 
generally applied to describe the behavior of conductance 
in a disordered wire. These two approaches give very 
similar solutions for the probability distribution of con- 
ductance in Q1D and predict some universal properties of 
electron transmission. They also give very similar expres- 
sions for the localization lengths (LL): £m m (M +1)1 and 
Cm ~ [/3{M — 1/2) + l]l, respectively (M is the number 
of the propagating channels, I is the phenomenological 
mean free path, and (3 — 1(2) for orthogonal (unitary) 
symmetry class). In either approach, however, the phe- 
nomenological I is viewed as a fixed parameter. The ques- 
tion how LL explicitly depends not only on energy E, but 
also on the coupling constants and the type of boundary 
conditions in Q1D disordered systems is left out in these 
analyses. What we are trying to point out is that in 
spite of the progress which has been made to towards a 
characterization of the localization in Q1D systems, mi- 
croscopic analytic studies of LL as a quantum parameter 
of fundamental importance, has still not been achieved. 

The first step in this direction was done by Dorokhov 
in Ref. Q, who calculated the LL of M random tight- 
binding (TB) chains with random site-energies. The LL 
in a weak disordered regime was obtained by the author 
for a Q1D wire with M channels, is independent of the 
number of channels A/. This result was questioned by 
Heinrichs [H, where it was shown that for weak disorder 
and for coupled two- and three-chain systems (M = 2, 3) 
the inverse LL is proportional to M, in contrast to the 
result of Ref. 4]. However, this approach, adopted in 
[H and based on a scattering matrix treatment of con- 
ductance, does not allow author to extend his studies of 
LL to Q1D systems with larger numbers of scattering 



channels M. Recently, progress has been made in tak- 
ing into account an arbitrary number of channels in the 
calculation of LL. Romer and Schulz-Baldes using a 
perturbative formula for the lowest Lyapunov exponent 
(the inverse LL) for Q1D TB Anderson strip model with 
PB conditions obtained LL's dependence on energy E, 
propagating modes M (M is even) and disorder strength 
W . A non-perturbative analytical approach, based on 
the Green's function formalism to solve the Dyson equa- 
tion in Q1D and two-dimensional (2D) disordered sys- 
tems without any restriction on the numbers of impu- 
rities and modes, was developed in Refs. F° r a 
TB Hamiltonian with several modes and on-site disorder 
the electron's scattering matrix elements T nm (hence the 
wire conductance G — J2 nm TnmT* m (in units of e 2 /h)) 
were analytically calculated for an arbitrary impurity 
profile without actually determining the eigenfunctions. 
In these papers only HW conditions were discussed, 
which correspond to arranging the parallel equidistant 
chains on a plane. 

We have performed a careful numerical analysis and 
have derived the LL in a Q1D system with HW and PB 
conditions. In our numerical calculations we used Kubos 
formula for computing conductivity. The numerical re- 
sults were compared with the existing analytical expres- 
sions of LL, calculated in Refs. |j-{9|. Surprisingly, our 
numerical calculations show that none of these expres- 
sions for LL fit the numerical data well. Particularly, 
LLs calculated in Refs. [3, H, 0, H| result in an incorrect 
dependence on M, while LLs calculated d, 0] correctly 
predicted the M dependence, but failed to provide the 
exact magnitude of LL. 

The discrepancy between the theory and our numerical 
results, as we will explain later on, is due to two main 
factors: (i) Since it is not easy to calculate the right- 
hand term of the equation (7) (see below), theoretical 
calculations assume that in the weak disordered regime 
the length (In G) can be replaced by length In (G) and 
by expanding to lowest order of the powers of the dis- 
order and, after averaging over realizations one can get 
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a closed analytical expression for LL in the Q1D system 
((...) denotes averaging over disorder realizations). How- 
ever, because of not self-averaging the conductance G, 
these lengths do not agree with each other and thus lead 
to a different answer for LL. Note that the same type of 
problem exists also in ID disordered system where, in the 
weak disordered limit, LLs, obtained numerically and an- 
alytically, differ by a factor of 2 (see e.g., Refs. [13, El), 
(ii) As follows from the numerical analysis of the rela- 
tionship between the different LLs in Q1D systems (see 
Eq. (1)), the right hand-side term is not zero. This is an 
essential piece of information, which allows us later on 
to introduce new LLs for different boundary conditions 
in the transverse direction, which fit the numerical data 
very well. 

It is worth noting that while in ID the relationships 
between the various lengths is well known (e.g., (InG) = 
41n(G) = -\ In (1/G)), in Q1D, to the best of our knowl- 
edge, no such calculations have been previously reported. 
Our first goal consists in checking numerically what rela- 
tionship exists for different lengths in the Q1D case. Once 
this is established, motivated by our doubts about the 
correctness of LLs results of Refs. and to overcome 

the difficulty of the discrepancy, we have reconsidered the 
calculation of LL for the Q1D TB anisotropic Anderson 
modeljUsing the Greens function approach, developed in 
Refs. 17-91. This is our second and main goal. The ana- 
lytical results for LL with HW and PB conditions, Eqs. 
(9) and (11), are then compared with numerical results. 
Excellent agreement with analytical calculations can be 
achieved if one multiplies the expressions (9) and (11) by 
a factor 2 and shifts them up by £i for HW and by £i/2 
for PB conditions, respectively. £i = 96sin 2 fci/W 2 is 
the LL in ID disordered system, calculated in the weak 
disordered regime [12j . 

Relationship between different lengths in Q1D sys- 
tems.— First we study the relationship between (InG) 
and In (G) in Q1D disordered systems. Our numerical 
calculations show, that these two lengths are connected 
through the relation 



(InG) -4 In (G) = G. 



(1) 



The constant G is different for HW and PB conditions 
and for each case is determined numerically. G tends 
to zero in the ID case as expected. The length L de- 
pendence of (InG) and 41n(G) is plotted in Fig. 1 for 
M = 10 and W = 1. The slopes of the two lines are the 
same within error bars. For each value of W we used L's 
that ensure that we are well inside the exponential decay 
(see Fig. 1). 

Some technical details follow: to obtain the mean val- 
ues (...) we have used 10 5 independent realizations of the 
disordered strip. Assuming a Gaussian form (as we have 
checked) for the logarithm of the conductance, z = In G, 
it is straightforward to show that the average of G is 
(G) = J M GP(G)dG = J M PQnG)dG. Solving the last 
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FIG. 1: Averaged logarithm of the conductance (InG) and 
logarithm of the average conductance In (G) as a function of 
the length L of the strip. 



integral and combining it with the numerical results for 
(In G) and its variance a 2 we get 



(G) = -cxp((lnG) + (i 2 /2)Erfc 



(InG) 



InM 



V2a 



(2) 

where Erfc is the complementary error function. 

Q1D TB anisotropic Anderson model.— Let us dis- 
cuss a Q1D disordered lattice of size L x M described by 
the standard TB anisotropic Hamiltonian with nearest- 
neighbor transfers, t x and t y along the x and y directions, 
respectively 

L M 

3 = 1 1=1 (3) 

-EE ; ) ** 0' + M + li, t y (j, i + 6\}, 

3,1 S=±l 

where \ j, I) is the atomic orbital at site (j, I) and ejj is the 
strength of the random potential at site (j, I), assuming 
it to be uniformly distributed in the interval (—3-, ^-)- 
The disordered region is connected to perfect leads on 
both ends, extended to ±00 in the x direction. L is the 
length of the system and M is the number of modes in 
the left and right leads. For simplicity we choose the lat- 
tice constant to be equal to 1. For further calculations 
we assume the existence of a confining potential V c (yi) in 
the discrete y direction (yi = I, I = 1, 2, M) leading to 
a set of transverse modes, whose actual values depend, 
however, on BC conditions. For HW and PB conditions, 
the energy of the electron is given by the following dis- 
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persion relations 



E = 



*2itj; cos 



Of rn s m% 



ti= 1,2,..., M, HW, 



2t x cosk n + 2t y cos^, n = 0,1, ...,M- 1, 



PB. 

(4) 

The appropriate eigenf unctions, ip n {yi), of the ID 
Schrodinger equation with periodic potential of the chain 
of atoms along the y direction with HW and PB condi- 
tions are (I = 1,2, ...,M) 



2_ ■ ■nln 
M Bm M+l 



HW, 



2-<5„ n-5 



M 



^exp[i^] 



(5) 



PB. 



Next, closely following Refs. we evaluate the 

scattering matrix elements T nm , in the weak disordered 
regime. The result for the electron transmission ampli- 
tude Tr, m IS 
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fe„ is the Fermi wave vector of the n— th subband (chan- 
nel) and is determined by the energy dispersion relation 
flU). For M = 1 it reduces to the LL £i for a ID chain. If 
there is no coupling to the second, third, etc, modes, all 
k n are equal, and after the summation over the modes, we 
find from Eq. ® £ M W = M£i. This result is somewhat 
expected: it confirms the prediction of Thouless [13] that 
in the limit of weak coupling £ M W must be proportional 
to M. Although one can get two correct limiting values 
£i and £m>i from expression ^ w , Eq. ([9]), it fails to 
give the exact value of LL for an arbitrary M. 

Our direct numerical computation of the LL for the 
Anderson model ^ shows that we can get an almost 
perfect agreement with the theoretical £ M W , Eq. © for 
M > 2, if we multiply the latter by a factor 2 and shift 
it up by £i, i.e., redefine new LL Xm W 



■v HW 

AM 



(HW 
SM 



M > 2. 



(10) 



where A t = i £ n=1 S '" s£fc„ ' ^ = ( fc « _ fc ™)^ _ X ) 
and D = 1 + i E/=i Ei=i The wave numbers /c n 

for the propagating modes are defined by Eq. (j4}, for HW 
and PB conditions, respectively. Similarly, L t is equal to 
M + l or M depending on the BC. 

The inverse LL £m as a function of the system size L 
and modes M can be written as 
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which is valid for both boundary conditions. 

Hard wall conditions.— Using the explicit expressions 
for ip n {yi) (see Eq. ©) and © for the LL £ M with HW 
conditions, when the M channels are propagating, we 



Now, replacing (In G) by In (G) and assuming that for 
weak disorder the transmission coefficients are close to 1 
and thus the reflection coefficients are close to zero, we 
can expand the right-hand side of Eq. ([7]). Next, after 
ensemble averaging over the random potentials eji dis- 
tributed uniformly according to the explicit expressions 
for T nm , Eq. ([6]) and keeping the terms to order W 2 , we i 
arrive at the following expression for the inverse LL Tpb 

SM 



Figure 2 shows the M mode dependence of \m ■ The 
solid lines have been computed from Eq. (TTU1) and the 
dots denote the result of numerical calculations. The 
good agreement found fully supports the validity of the 
analytical expression for Xm • The comparison is free 
of any adjustable parameter. We have checked that the 
analytical expression, Eq. (fT0|) , agrees very well with the 
numerical data in the whole range of vertical hopping 
parameter < t y < 1 where the Xm ^ s a linear function 
with respect to M. The non-linearity starts when t y > 1 
and for those t y the numerical and analytical results start 
behaving differently due to the fact that the validity of 
formula (fTUl) breaks down. For these t y > 1 our data 
show that L/xm W ~ 1; i- e - takes place at the crossover 
from a Q1D to 2D system. 

Periodic boundary conditions.— The result for the LL 
reads 



96 



i 

M 2 



M 



2 — &i,o~ S 



(M-l) 3 



i=0 sin ki 
M- l 

E;=o 



if M even 



(11) 



2-5, 



'1= sin ki 



if M odd, 



where k n must be defined from the dispersion relation 
©. 

The process of deriving the expression for even M is 
quite straightforward. Using the explicit form of the elec- 
tron wave function ip n (yi), (0) and Eq. jS} yields the 
desired result. The case for odd M requires special con- 
sideration. First, for the infinitely long periodic system 
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FIG. 2: The dependence of the localization length Xm on 
the number of modes M for disorder W = 0.4 and 1.0. Dots 
are the numerical results and each data point corresponds to 
an average over 10 s realizations of disorder. The solid lines 
represents the theoretical prediction, Eq. (|12|l . At E — 0.02 
all modes are propagating. 



the conductance G is an odd function of energy, which 
is in contrast to the symmetric behavior of G with even 
M modes. Second, the analysis of the conductance of 
the ideal TB model as a function of the energy (at fixed 
even M) shows that the change from one plateau value 
to the next one is 2 (in units of e 2 /h), while in the case 
of even M, it is 1. Formally this means that M must be 
replaced by (M — 1) in the expression of LL with even 
number of M . This conjecture was numerically tested 
and supported by the direct numerical calculation of the 
LL (see Fig. 3). It is clear that the difference between 
M and (M — 1) is negligible for large M, but may not be 
negligible for small M . 

As in the case of HW condition we get an excellent 
agreement with the theoretical £,^j B , Eq. (fTTj) for M > 2, 
if we multiply the latter by a factor 2 and shift it up by 
£i/2. The new LL xli is 
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pPB 
SM 
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M > 2. 



(12) 



In Fig. 3 we have tested the prediction of the analytical 
theory against the numerical results where the M mode 
dependence of x™, Eq. (fl~2"j) is shown. Solid lines have 
been computed from Eq. (fT2")l and dots denote result 
of numerical calculations. The good agreement between 
simulations (dots) with Eq. (fT2"j) is evident for a relatively 
large range of disorder W. In the right panel of Fig. 3 
our numerical data for LL was compared with similar ex- 

(dashed line) from 



w 

pression | ~ ^ l m=Q sin kl sin kji 
Ref. [6|- One can see that the slope of the dashed line 
agrees with numerical calculations, but certainly there is 
a problem with accurate numerical values of LL. To get 
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FIG. 3: M dependence of the localization length defined by 
Eq. (|12[) . Dots are the numerical results. Left: M is odd. 
Right: M is even. Dashed line is presents Eq. (8) of Ref. 



a correct value for LL for an arbitrary M one needs the 
dashed line to shift up by about 0.39£i. 

In summary, our numerical calculations indicate that 
in Q1D disordered systems the right-hand side term of 
Eq. ([1]) is not zero. This means that without a vertical 
shift it is impossible to explain the discrepancy between 
the theory and numerical data and get the correct an- 
alytical expression for LL in the framework of different 
approaches used in Refs. Based on Eq. ([TJ we 

have presented the analytical expressions (fT0|) and (fT2|) 
for LL which are in excellent agreement with numerical 
calculations. 
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